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THERMODYNAMIC PROPERTIES OF CESIUM UP TO 1500° K 

by Sheldon Heimel 
Lewis Research Center 

SUMMARY 

Consistent tables of thermodynamic properties of elemental cesium were compiled 
for both the pure species and the equilibrium vapor on the saturation line using selected 
values of 10 500 calories per mole for the heat of dissociation of dimer and -18 920 cal- 
ories per mole for the heat of condensation of monomer. The equilibrium vapor prop- 
erties and properties of condensed cesium (Cs) are given to 1500° K, while properties of 
gaseous monomer and dimer are given to 2500° K. 

The table of mixture properties on the saturation line gives the following thermo- 
dynamic functions: temperature T, total vapor pressure P, weight fraction of dimer in 
vapor X 2 , molecular weight Mg, specific volume v, density p, enthalpy of condensed 
phase he, enthalpy of vaporization Ah^p, enthalpy of vapor h g , entropy of condensed 
phase s c , entropy of vaporization As ya p, entropy of vapor Sg, and heat capacities 
(equilibrium and frozen) (c p ) eq and (c p ) fr . Assigned enthalpies are relative to enthalpy 

of the crystal at 298. 15° K, which is taken to be zero. The tables of properties for the 
pure species give the following thermodynamic functions for the standard state: heat 
capacity at constant pressure C°, sensible enthalpy H^, - H^, entropy S^,, sensible 
free energy F^, - Hq, the sum of sensible enthalpy at T° K and chemical energy at 
0° K H^,, and values of enthalpy changes and logarithms of the equilibrium constants. 
The latter two functions are given for the reactions of formation of Cs(c), Cs^(g), and 
Cs 2 (g) from the assigned reference state Cs(c) ((AH^) and log^Q Kj) and from Cs^(g) 

(AH^, and log 1Q K). 

Thermodynamic functions for the gases were generated from atomic and molecular 
data, whereas the functions of the condensed phase were based on selected experimental 
data that were smoothed and made self-consistent. 


INTRODUCTION 


Current work at the Lewis Research Center indicates the need for reliable vapor 
pressure data as well as other thermodynamic properties for cesium in the following 
areas: 

(1) Rankine cycle engine. The thermodynamic properties of cesium are needed in 
order to make an analysis of a Rankine cycle system and its components. 

(2) Thermionic plasma diode. Improved performance requires knowledge of the con 
centration of cesium in the gas phase and the rate of evaporation from a liquid surface. 
Such data can be derived with the help of an equation for pressure as a function of satura 
tion temperature. 

Accordingly, to meet these needs, the following tables have been generated: 

O 

(1) Boiling points corresponding to pressures of from 10” to 10 atmospheres 

(2) Properties on the saturation line from 100° to 1500° K: temperature T, °K; 

total vapor pressure P, atmospheres; partial pressure of monomer p^, 
atmospheres; partial pressure of dimer p 2 , atmospheres; weight fraction of 
dimer x 2 , gram per gram mixture; molecular weight of vapor Mg, gram per 
mole; specific volume of vapor v, cubic centimeter per gram; density of 
vapor p, gram per 1000 cubic centimeters; enthalpy of condensed phase h c , 
calories per gram; enthalpy of vaporization Ah^p, calories per gram; en- 
thalpy of vapor h , calories per gram; entropy of condensed phase s , cal- 

& Q c 

ories per gram per K; entropy of vaporization As ya p, calories per gram per 
°K; entropy of vapor Sg, calories per gram per °K; equilibrium heat capacity 
of vapor (c ) , calories per gram per °K; and frozen heat capacity of vapor 

p eq „ 

(c ) , calories per gram per K 

p fr 

(3) Thermodynamic functions for each pure species: temperature T, K; heat 

capacity at constant pressure C°, calories per mole per °K; sensible enthalpy 

H T " H 0’ ca l° r * es P er mole; entropy S^,, calories per mole per °K; sensible 
free energy F^, - H^, calories per mole; sum of sensible enthalpy at T° K 
and chemical energy at 0° K H^,, calories per mole; and values of enthalpy 
changes and logarithms of the equilibrium constants. The latter two functions 
are each given for two different reactions of formation. One set of functions, 
(AH^) and logjQ K f is for the reaction of formation of Cs^g) and Cs 2 (g) 


from the assigned reference state Cs(c). The other set of functions, AH^ 
and log 10 K, is for the formation of Cs(c) and Cs 2 (g) from Cs^g). The sym- 
bols g and c refer to the gas and condensed phases (crystal and liquid), 
respectively. The thermodynamic functions for the gases were generated from 
atomic and molecular data, whereas the functions of the condensed phase were 
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based on selected experimental data that were smoothed and made self- 

consistent. The functions for the condensed phase are given up to 1500° K; 

those for the monomeric and dimeric gases are given up to 2500° K. 

When vapor pressure data are available, they may be combined with sensible free 

energy data by use of the third law of thermodynamics to calculate an implied value of 

(AH^) (heat of condensation of monomer at 0° K) for each available vapor pressure 
c 

point. If the vapor pressure data are completely consistent with the thermodynamic 
data, the calculated values of (AH^) will be constant. The usual result, however, is 

that there is variation in the calculated values of (AH^J . Depending on the extent and 

u c 

nature of this variation, some average (AHq) may be selected. This selected value of 

u c 

(AHq) together with the sensible free energy data may now be used to calculate vapor 
u c 

pressures. These calculated vapor pressures and the experimental vapor pressures 
will, of course, not agree. At this point either the experimental vapor pressures must 
be abandoned or the theoretical free energies must be modified. 

In this report a similar procedure to that just outlined will be followed. There are 
several differences in the analysis, namely, 

(1) The experimental vapor pressures are first smoothed. 

(2) The gas is assumed to be a mixture of monomer and dimer. Therefore, the 
intermediate step of calculating equilibrium compositions is required. 

(3) The procedure for selecting a value of ( AHq) is to find that value of ( AHq) 

u c ^ c 

which gives a minimum deviation of calculated and smoothed experimental vapor pres- 
sures over the entire range. 

(4) The selected (AH^.) and thermodynamic data are used to calculate a recom- 

u c 

mended set of vapor pressures for the crystal and the liquid. 


PRESENTATION OF TABLES 

Table I contains the smoothed values of C° for crystalline cesium up to 100° K 

r 

together with consistent values of H° - Hq, S^,, and -(F^ - Hq). Generation of these 
data will be discussed in SELECTION OF INITIAL DATA. 

Tables n to IV (pp. 5-7) give the thermodynamic properties of Cs (crystal, liquid) 
up to 1500° K; Csj(g) up to 2500° K; and Cs 2 (g) up to 2500° K, respectively. The 
properties presented are C°, H° - h£, S° , -(F° - H°), H° , (AH°) f , log 1Q Kj, AH° , 

and logjQ K. In table II, data for the crystal extend up to the melting point, 301. 8° K. 
Data for the liquid are tabulated at 100° intervals from the melting point to 1500° K. In 
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TABLE I. - THERMODYNAMIC PROPERTIES 


OF CESIUM (CRYSTAL) UP TO 100° K 


T, 

°K 

C° 

cal/(mole)(°K) 

TT° TjO 

cal/ mole 

S° 

cal/(mole)(°K) 

-(F° - H°), 
cal/mole 

5 

0.746 

1.2 

0.3342 

0.5 

10 

2.379 

9.0 

1.3539 

4.5 

15 

3.719 

23.9 

2.5835 

14.9 

20 

4.680 

44.5 

3. 7993 

31.5 

25 

5. 110 

69. 1 

4.8935 

53. 3 

30 

5.385 

95.3 

5.8513 

80. 2 

35 

5.558 

122.7 

6.6955 

111.6 

40 

5.668 

150.8 

7.4454 

147.0 

45 

5.741 

179.4 

8. 1174 

185.9 

50 

5.795 

208. 2 

8. 7252 

228. 1 

55 

5.842 

237. 3 

9. 2798 

273. 1 

60 

5.888 

266.6 

9. 7901 

320. 8 

65 

5.933 

296.2 

10. 2631 

370.9 

70 

5.977 

325.9 

10. 7044 

423.4 

75 

6.016 

355.9 

11. 1182 

477.9 

80 

6.049 

386. 1 

11.5076 

534. 5 

85 

6.075 

416.4 

11.8751 

593.0 

90 

6.097 

446.8 

12. 2230 

653. 2 

95 

6.122 

477.4 

12. 5533 

715. 2 

100 

6.160 

508. 1 

12.8682 

778. 7 


tables HI (Csj) and IV (Cs 2 ), data are 
tabulated at 100° intervals from 0° to 
2500° K for the ideal gases. Table V 
(p. 8) contains the molecular con- 
stants used to calculate thermody- 
namic functions for dimeric cesium. 

Table VI (p. 8) contains recom- 
mended boiling points of liquid 
cesium calculated from the final set 
of thermodynamic properties. These 

boiling points correspond to pres- 
-8 

sures from 10 to 10 atmospheres 
at every power of 10. For compar- 
ison, boiling points obtained from 
equation (33) are also given. 

Table VII (p. 9) gives the follow- 
ing thermodynamic properties of the 
equilibrium gaseous mixture on the 
saturation line: T, P, pp p 2 , x 2 , 

Mg’ P. h c> ih vap’ h g’ s c’ as vap' 

s , (c ) , and (c ) . The values of 

g p eq p fr 

P, Pp P 2 > v > anc * P are given in 
floating-point notation, where the 
decimal number is to be multiplied by 


10 raised to the power of the sign and the two digits following the letter E; for example, 
0. 100000E-03 is 0. lOOOOOxlO -3 = 0. 0001. 


Assigned Reference State 

The assigned reference state is Cs (crystal, liquid). The crystal is the reference 
state up to the melting point, and the liquid phase is the reference state above the melt- 
ing point. 


Assigned Enthalpy Values 

For some applications (ref. 1), it is convenient to combine sensible enthalpy 
H^, - Hq and chemical energy into one numerical value, H^. The arbitrary base 

4 




TABLE II. - THERMODYNAMIC PROPERTIES OF CESIUM (CRYSTAL, LIQUID) 
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TABLE m. - THERMODYNAMIC PROPERTIES OF Cs 1 (GAS) 
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The reference element is crystal cesium up to 301. 8° K and liquid cesium above 301. 8 ( 









TABLE IV. - THERMODYNAMIC PROPERTIES OF Cs 2 (GAS) 
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The reference element is crystal cesium up to 301. 8° K and liquid cesium above 301. 8 ( 












TABLE V. - MOLECULAR 


TABLE VI. - BOILING POINTS OF LIQUID 


CONSTANTS FOR DI- 


a Constant from ref. 10. 

^Constant from ref. 11. a log 10 P atm" ~ 4 °y — " + 7 - 04453 - °- 915282 log 10 T. 


ATOMIC CESIUM 


Molecular weight 

265.82 

Symmetry number 

2 

Electronic state a 


Statistical weight 

1 

-p 

w e’ cm a 

41.990 

o) x , cm" 1 

e e - i a 

^„y p , cm 

_ e _ib 

B , cm 1 

0.08005 

-0.0001643 

0.01272 

-l b 

a , cm 1 

0.000035 

6 i b 

D , cm" 1 
e 

0.467X10" 8 


CESIUM FROM 10‘ 8 TO 10 ATMOSPHERES 


lo Sl0 P atm 

Boiling points, °K 

Recommended 

From equation (33) a 

-8 

317. 8 

317.8 

-7 

345.8 

345.8 

-6 

379. 2 

379.4 

-5 

420. 1 

420. 3 

-4 

471. 2 

471.4 

-3 

537. 0 

537.2 

-2 

625.0 

625.0 

-1 

749. 1 

748.7 

0 

937. 7 

937. 3 

1 

1264.0 

1264.5 


for assigning values to the enthalpy of Cs(c) was a value of zero at 298. 15° K. Since 

table n gives (H^ 9g 15 - H^) = 1807. 2 calories per mole, then (Hq) = 

Cs(c) Cs(c) 

-1807. 2 calories per mole. From table II, (AhS.) = (Hji) - (h2.) = 

° c ° Cs(c) ^ C Sl (g) 

-18 920.0 calories per mole. Therefore, (h2.) = -1807. 2 - (-18 920.0) = 

° c Sl (g) 

17 112. 8 calories per mole. From table IV, (AH2.) = -D2. = (h2.) - 2(h£.) 

° 2 OX) Cs2 ( g ) X) CSi ( g ) 

-10 500 calories per mole. Therefore, 


(Ha) = -10 500. 0 + 2(17 112. 8) = 23 725. 6 cal/mole 

^ Cs 2 (g) 


Heats of Formation 

Two sets of values for heats of formation are given in tables II to IV. The first set 

(col. 7) is for the formation of the given species from Cs(c), (AIlS,) ; the second set 

1 f 

(col. 9) is for the formation from CSj(g), (AH^,). For Cs 2 (g) at 298. 15° K, for ex- 
ample, 


( AH 298.15) - ~ ( H ° 


f, 2 298 ‘ 15) Cs 2 (g) " 2(H298 - 15) Cs(c) 


U) 
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TABLE VH. - THERMODYNAMIC PROPERTIES OF CESIUM ON SATURATION LINE 
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( AH 298 15> = < H 298 15 > “ 2 ( H 298 15 ) 

298. ID 2 298.15 Cg2(g) 298. 15 Cg (g) 


( 2 ) 


Equilibrium Constants 


Two sets of logarithms of the equilibrium constants for the two formation reactions 

discussed previously are also listed in tables II to IV. The equilibrium constant Kj for 

formation from the assigned reference element (col. 8) is obtained from the standard 

free energy change ( AF™) by means of the equation 

1 f 


(AF T } f 

l0g l° ^ 2. 3025851 RT 


(3) 


The equilibrium constant K for formation from the atomic gas is obtained from a sim- 
ilar equation 


log 10 K= -- 


AF" 

2. 3025851 RT 


(4) 


Atomic weights, the universal gas constant and the constants used in the evaluation 
of the entropy constant were the same as those used in reference 1. 


SELECTION OF INITIAL DATA 
Crystal 

Heat capacities of the crystal (table I, p. 4) were derived by smoothing the experi- 
mental data from the sources given in table VIII. 

The following procedure was used to obtain smooth data from 0° to 100° K. A 
least-squares fit of the data of reference 2 was used in the temperature range from 0° to 
4° K. Data from 4° to 20° K were obtained from a curve drawn through the data of ref- 
erence 2 at 4° K and through the data of references 3 and 4 at 20° K. Data from 20° to 
100° K were obtained from a least-squares fit of the data of reference 4 in that tempera- 
ture range. 

Between 100° and 210° K, the data of reference 4 are approximately linear in tem- 
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TABLE Vm. - REFERENCES FOR HEAT 
CAPACITIES OF CRYSTAL 


perature. The heat capacities at 100° and 
210° K yield the equation 


Temperature 

Source 

range, 

°K 


0. 3 to 4 

Ref. 2 

8 to 12 

Private communication from 


Professor D. C. McCollum 
of University of California 
in Riverside, California 
and ref, 3 

20 to 210 

Ref. 4 


C°= 5. 76 + 0.004 T (5) 

Above 210° K, thermal effects associated with 
the cesium-oxygen system cause anomalous in- 
creases in heat capacity and therefore these 
data were not used. Equation (5) was used to 
extrapolate heat capacities up to the melting 
point. The melting point of 301. 8° K and the 


heat of fusion of 520. 1 calories per mole were taken from reference 5. 

Reference 4 gives a melting point of 300. 5° K and the total heat absorption between 
the crystal at 262. 5° K and the liquid at 300. 5° K. If this heat is compared with the 


sensible enthalpy change of the crystal between these two temperatures, as calculated 
from equation (5), the implied heat of melting is 514 calories per mole. This value is 


in reasonable agreement with the value adopted. 

The enthalpy from the linear equation for C°, when combined with the heat of fusion, 

r 

gave 




( 6 ) 


Liquid 

A number of empirical equations have been used to satisfactorily fit enthalpy data as 
a function of temperature for various substances; however, heat capacity data derived 
from differentiation of the enthalpy equations are often very unreliable. For any par- 
ticular species, various empirical enthalpy equations that might be chosen yield con- 
siderably different sets of heat capacities. 

Two sets of experimental enthalpy data were considered: that of references 6 and 7. 
The curve of reference 6 exhibits an anomalous bump from 351° to 620° K, followed by 
a straight line to 1238° K. It is recommended in reference 6 that just the straight-line 
portion be used, which has a standard deviation a = 1. 5 percent. For the cubic equa- 
tion (from 340° to 1176° K) of reference 7, a = 1. 1 percent. 

The implied C° = 7. 25 calories per mole per °K of reference 6 differs greatly 
from the heat capacities derived from the enthalpy curve of reference 7 and listed in 
table IX. 

According to reference 7, the curves of heat capacity against temperature of liquid 
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TABLE IX. - HEAT CAPACITIES 
DERIVED FROM ENTHALPY 
CURVE OF REFERENCE 7 


Absolute 

temperature, 

T, 

°K 

Heat capacity at 
constant pressure for 
standard state, 

<$• 

cal/(mole)(°K) 

301.8 

10.35 

500 

7. 28 

700 

6. 23 

1000 

8. 53 

1200 

12.64 

1500 

22. 67 


cesium and other alkali metals have a parabolic shape; 
that is, heat capacity initially decreases with increas- 
ing temperature, reaches a minimum, and then in- 
creases continuously. However, a parabolic shape for 
heat capacity of these species is only the result of 
selecting a cubic to represent experimental enthalpy 
data. 

It was thought advisable to use a straight line for 
enthalpy because (1) the straight line, in the standard 
deviation sense, represents cesium enthalpies about 
as well as the cubic does and (2) the cubic has the dis- 
advantage of implying very large heat capacities for 
even a few hundred degrees extrapolation. 

The actual straight line used in this report was 
based on the straight-line data of reference 6 but was 


constrained to go through the previously calculated relative heat content of 545. 5 cal- 


ories per mole at the melting point (eq. (6)). This gives the following equation: 


(H°) 

1 Cs( liquid) 


“ (Hoqq 15) 

^ yy - lb Cs(crystal) 


= -1704.5 + 7.455 T 


(7) 


Equation (7) implies a constant value of C° = 7. 455 calories per mole per °K. 

Had the data of reference 7 been fitted to a straight line and similarly constrained to 
go through the heat content of 545. 5 calories per mole, a constant value of 0° = 7. 595 
calories per mole per °K would have resulted. The enthalpies would be, on the average, 
only about 1 percent higher than those derived from the data of reference 6. 

The enthalpy data of reference 8 are about 6 to 7 percent lower than the data of both 
references 6 and 7 over the range from about 600° to 1300° K; therefore, they were not 
used. 


Monomer (Gas) 

The partition function used to compute the thermodynamic properties of Cs^ was 
truncated by the temperature -dependent cutoff technique used in reference 1. The re- 
sults of these computations are exactly the same up to 2200° K as would be obtained by 
simply using all levels given in reference 9 with no cutoff. Above 2200° K enthalpy and 
entropy differ only in the last place computed and heat capacity in the last or second last 
place. 
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Dimer (Gas) 


The method used to compute the thermodynamic properties of CS 2 is the method 
for diatomic molecules described in reference 1. The constants for the dimer are 
shown in table V (p. 8). They were taken from references 10 and 11. 


Vapor Pressures 

During the period 1913 to 1937, the investigations reported in references 12 to 18 
obtained vapor pressure data for cesium in the temperature range from 238° to 670° K. 
This early work was evaluated in reference 19 where it was concluded that, except for 
the results of references 17 and 18, the early work is not too reliable. 

The vapor pressures in the 455° to 589° K range (ref. 17) were used in reference 19 
to generate a vapor pressure equation. Inasmuch as the data of reference 17 are only 
relative pressure data (obtained by the magneto -optical method) and are converted to 
absolute pressure data by relating them to the unreliable data of reference 13, it seems 
prudent to disregard the data of reference 17. 

Recently, the first experimental measurements in the temperature range that in- 
cluded the boiling point (742° to 1199° K) were reported in reference 20. These data 
check reasonably well with the later data of reference 6 (727° to 1334° K). 

In view of the previous observations, the data of references 18 (vapor pressures in 
the liquid range), 20, and 6 were combined and fitted by the least-squares technique. 

All the data were given the same weight. The results of several empirical equations 
were compared. The following equation, which is in the form of the Kirchhoff equation, 
was accepted because it gave as satisfactory results as the other forms considered: 

log 10 P = - 392 °- 3 - + 5. 71342 - 0. 519781 log 10 T (8) 

If only vapor pressures were considered in this report, then the smoothed experi- 
mental vapor pressures obtained from equation (8) would be accepted. As pointed out in 
the INTRODUCTION, however, these smoothed data, when combined with other thermo- 
dynamic data, lead to inconsistencies in the heat of condensation of the monomer at 0° K 

(AH^) . Therefore, the data from equation (8) are used only for preliminary analysis, 
c 

whereas the final recommended vapor pressures are calculated as explained in detail in 
the section Pressures and Weight Fractions of Gaseous Species. 
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Heat of Dissociation at 0° K 


An upper bound (4020 cm" *) for the heat of dissociation at 0° K(Dq) is estimated in 
reference 21 by analytical extrapolation of the ground state vibrational intervals. A 
lower bound (3318 cm”'*') is estimated in reference 21 by making certain assumptions as to 
the final states of the excited atoms produced. The average of these two values, which is 
equivalent to 10 500±1000 calories per mole of dimer, is recommended in this reference. 

The estimate of 0. 453 eV or 10 450 calories per mole from the correlation of force 
constants, equilibrium internuclear distances, and heats of dissociation of the alkali 
metals is obtained in reference 22. A technique for deducing a value of Dq where Dq 
was not well known is suggested in reference 23. This technique gave an optimum 
Dq = 10 000 calories per mole which is in fairly close agreement with the spectroscopic 
value of 10 500 calories per mole. 

The calculation of (AHq) , which is discussed in the next section, depends in part 

c 

on the value of Dq. While no value of Dq that was tried gave a constant value of 
(AHq) over the entire temperature range of 301. 8° to 1500° K, values of Dq from 

about 10 000 to 11 500 calories per mole gave the most nearly constant value of (AH^) . 

Based on this analysis, the conclusion was that the spectroscopic value was ap- c 
proximately correct, and the value Dq = 10 500 calories per mole was selected. 

Heat of Condensation of Monomer at 0° K 

The values of (AH2.) in calories per mole are given in a number of publications as 
w c 

follows: -19 048 (ref. 24), -18 790 (ref. 19), -19 050 (ref. 25), and -19 035 (ref. 26). 

The aforementioned values of (AH^.) are based on the pre-1962 vapor pressures. 

c 

Therefore, this report will recommend a "best" (AHq) of -18 920 on the basis of the 
newer vapor pressures, as discussed in the following section. 


THERMODYNAMIC ANALYSIS 

As discussed in the INTRODUCTION, the procedure for the selection of a value for 

(AH^) is to find that (AH^) for which the difference between the calculated and 
c c 

smoothed experimental vapor pressures is a minimum. As a first step, it is necessary 

to compute (AH$.) for each temperature in the range of interest. 

^ c 
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Derivation of (AH$) 


c 


The zero-degree heat of condensation of monomer was computed by means of the 
relation 


(AH^) = (AF°) - A(F° - Hq) 


where 


(AF°) = (F°) -(F°) =RTln Pl 

T c T Cs(c) r CSj(g) 


and 


A(F° - H^) = (F° - id) - (F° - id) 

T ^c T °Cs(c) T °C Sl (g) 


( 9 ) 


( 10 ) 


(ID 


Partial Pressure of Monomer 

The partial pressure of monomer p^ in atmospheres is obtained from the simul- 
taneous solution of Dalton’s Law and the equilibrium -constant equation for the dimeriza- 
tion reaction 


2 Csj = Cs 2 


The equations are 


P - + p 2 


( 12 ) 


and 



(13) 


where P is the total pressure, p 2 is the partial pressure of dimer , and K 2 is the 
equilibrium constant of dimerization. Values of K 2 are obtained from table IV (p. 7). 
The solution of equations (12) and (13) is 


15 


(14) 


?1 = 


-1 + V4PK 2+ 1 

" 


Therefore, apparent values of (AH~) can be computed at every temperature of interest 

u c 

by starting with given values of Dq, values of F^, - for monomer, dimer, and 
condensed phase, and smoothed experimental vapor pressures P F j T . This procedure 
can be summarized as follows: 


(F° - fig) 


Cs(c) 
(Table II) 


(F° - Hq) 


Cs x (g) 
(Table IH) 


Equation (11) 


A(F° - ig) 


K 2 

(Table IV) 


Equation (9) 


(Alg) 


FIT 

(Eq. (8)) 


Equation ( 14) 


Pi 


RESULTS OF CALCULATING (AH&) 

u c 

The apparent (Alg) values computed by means of this procedure from equation (9) 
u c 

are given in table X. The maximum difference is about 1 percent. 

If some single arbitrary value of (AH2.) is selected, then for each of the previous 

u c 

temperatures there will be a difference between the smoothed experimental vapor pres- 
sure Ppjr and the vapor pressure Pq^lC calculated with this arbitrary (AHq) 
value. c 

For this analysis, an error function may be defined as | P Frr - PcALc|/ p FIT 
averaged over the previous number of temperatures. The procedure for selecting a value 
of (AHq) c is to find that value of ( AHq) c which gives the minimum error function. 

The results of calculating the error function for various values of (AH&) c are given 
in table XI. The optimum (AHq) c was taken to be -18 920 calories per mole to the 
nearest 5 calories. 

The fact that there is still about 4 percent error between the smoothed experimental 
vapor pressure and vapor pressure calculated with the optimum (AHq) c> indicates that 
16 



TABLE X. - APPARENT VALUES OF -(AH^ 
AS A FUNCTION OF TEMPERATURE 


Tempera- 

ture, 

T, 

°K 

*( AH o> c > 

cal/mole 

Tempera- 

ture, 

T, 

°K 

-< AH oV 

cal/mole 

301.8 

18 791. 6 

1000 

18 965. 1 

400 

18 872. 5 

1100 

18 954. 3 

500 

18 921.7 

1200 

18 937. 0 

600 

18 949. 3 

1300 

18 913.1 

700 

18 963.6 

1400 

18 881.9 

800 

18 969. 7 

1500 

18 843.4 

900 

18 970.0 




TABLE XI. - ERROR FUNCTION 
FOR VARIOUS VALUES 
OF -(AHq) 


-(AH£) c , 

cal/mole 

Error function, 

|ap|/p 

18 850 

0.0538 

18 900 

.0397 

18 910 

.0381 

18 920 

.0370 

18 925 

.0370 

18 930 

.0372 

18 950 

.0389 

18 970 

.0458 


discrepancies still remain which may be in any of the following data: 

(a) Experimental vapor pressures 

(b) D° 

(c) Molecular constants of dimer (and, consequently, its thermodynamic properties) 

(d) Heat content (and other derived thermodynamic properties) of the condensed phase 

(e) (AH°) c 

The 4 percent error may also be due, in part, to the presence of a higher polymer, for 
example, trimer or tetramer; however, this possibility has not been considered in this 
report. 


CALCULATION OF EQUILIBRIUM MIXTURE PROPERTIES 
ON SATURATION LINE 

As explained in the INTRODUCTION, selection of an optimum value of (AHq) makes 

it possible to compute a set of vapor pressures consistent with the sensible free energies 
of the pertinent species. These vapor pressures and other equilibrium properties on the 
saturation line (e. g. , enthalpy, specific heat, and entropy) are given in table VII (p. 9). 
The equations used to obtain these properties are given in the following sections. 


Pressures and Weight Fractions of Gaseous Species 


By starting with the optimized (AH2.) and the sensible free energies of monomer 

u c 

and condensed phase, monomer pressure was computed with equations (9) to (11). 
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was 


Dimer pressure was computed from equation (13). The total pressure P^alC 
obtained from the sum of p^ and P 2 » 

The weight fraction of monomer x^ was obtained from the mole fraction of 
monomer Np thus, 


N l = ~ 
1 P 


( 15 ) 


N, 


X 1 “ 


2 - N, 


(16) 


The weight fraction of the dimer X 2 was obtained from the defining equation 


Xf + x 2 = ! 


(17) 


Enthalpy of Vaporization and Entropy of Vaporization 

The heat of vaporization in calories per gram from the condensed state to the 
equilibrium mixture on the saturation line was obtained from 

_ 2 Xl (AH^) { x + x 2 (AH^) f 2 

ih -e = ^ 

where M 2 is the molecular weight of dimer (265. 82 g/mole). The heats of formation 
are obtained from tables III and IV (pp. 6 and 7). 

The entropy of vaporization in calories per gram is obtained from equation (18) by 
means of 


As 


vap 


Ah vap 

T 


(19) 


Enthalpy, Specific Heat, and Entropy of Cesium Vapor 


The enthalpy 


‘g 


in calories per gram of an equilibrium mixture can be defined by 
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( 20 ) 



i=l i=l 


Equation (20) for saturation conditions gives the same results as the usual equation for 
the enthalpy of a mixture 


h„ = h, + h 

g vap c 

where h is the enthalpy of the condensed phase in calories per gram (table VII, p. 9). 

The equilibrium specific heat at constant pressure in calories per gram per °K is 
defined by 



Using equation (20) results in equation (21) becoming 



or 




+ < c p) 


r 


From equation (17), it follows that 



Therefore, 




(h 2 - hj) - 


^*2 

Ot> 


(ah2,) 
r 2 

m 2 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 
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Values of (A H™) are tabulated in table IV (p. 7). Values of (3x 0 /3TL maybe 
1 2 * P 

calculated from the following equation: 


(AH T ) 2 x 2 (2 - x 2 )(l - x 2 ) 
RT 2 2 


Equation (26) may be derived from the following expression, which is equivalent to 
equation (13): 

x o ( 2 - x 2 ) 

rjr 6 & i c\r 


4P(1 - x 2 )‘ 


and from the van't Hoff isobar, namely, 


3 In K, 


(AH°) 


p RT 


The entropy of the gas mixture may be obtained from 


s g = As vap + s c 


where s c is the entropy of the condensed phase in calories per gram per °K (table VII, 
p. 9) and As va ^ is obtained from equation (19). Equation (29) gives the identical 
results as the usual equation for the entropy of a mixture. 


Molecular Weight of Mixture 

An expression for the molecular weight of the mixture M e (in g/ mole) in terms of 

o 

the molecular weight and weight fraction of the dimer may be obtained as follows: 


M = - = = * 

n n 1 + n 2 

M l + M 2 


2 - Xr 
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where n, n^, and ng are the number of moles of mixture, monomer, and dimer per 
gram of mixture, respectively. 


Specific Volume and Density of Mixture 

The specific volume v in cubic centimeters per gram may be obtained from the 
ideal gas law 


_ RT 82. 05971 T 

PM„ PM_ 

g g 

The density p in grams per 1000 cubic centimeters is, therefore, 

„ 1000 
p = 

v 


(31) 


(32) 


DISCUSSION OF VAPOR PRESSURES 
Vapor Pressures Above Liquid 

-8 

Table VI contains boiling points of liquid cesium at saturation pressures from 10” 
to 10 atmospheres at every power of 10. These data were calculated in two ways: 

(1) From the recommended thermodynamic data, as explained in the section 

THERMODYNAMIC ANALYSIS 

(2) By a least-squares equation fitted to the recommended P^LC va * ues ^h a 

standard deviation of 0. 5 percent 

logiQ P atm = + 7. 04453 - 0. 915282 log 10 T (33) 

On the average the two sets of boiling points differ by about 0. 2° K. It should be 
emphasized that the vapor pressure equation (eq. (33)) is recommended rather than 
equation (8). 

Figures 1 and 2 compare the experimental vapor pressures of liquid cesium with a 
curve drawn through the recommended vapor pressures. The agreement of the two sets 
is seen to be very good. Figure 1 covers the temperature range from 303° to 1334° K. 
Figure 2 covers only the range from 729° to 1334° K to facilitate detailed comparison. 


21 






Figure 3. - Comparison of experimental and calculated vapor pressures above crystalline cesium. 


Subsequent to the completion of the analysis, a vapor pressure equation in the tem- 
perature range 928° to 1558° K was presented in reference 27. On the average, these 
vapor pressures are about 2 percent lower than those recommended in this report. 


Vapor Pressures Above Crystal 

Experimental vapor pressures above the crystal are available from reference 18 at 
temperatures from 238° to 298° K. In order to be consistent with the thermodynamic 
data adopted for this report, however, vapor pressures from 100° to 301. 8° K were cal- 
culated in the same manner as were the vapor pressures above the liquid. 

On the average, the calculated vapor pressures between 238° and 298° K were 
30 percent lower than the experimental ones (fig. 3). This difference is not significant, 
however, inasmuch as at the low pressure levels involved, the average difference is 
only in the order of 10”^ atmosphere. Furthermore, the lower calculated vapor pres- 
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sures are consistent with the assumption of reference 26 that the experimental vapor 
pressures are too high. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 26, 1965. 
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APPENDIX - SYMBOLS 






rotational constant, cm" 1 

heat capacity at constant 
pressure for standard 
state, cal/(mole)(°K) 

equilibrium specific heat at 
constant pressure of mix- 
ture, cal/g 

frozen specific heat at con- 
stant pressure of mixture, 
cal/g 

reactive specific heat at con- 
stant pressure of mixture, 
cal/g 



TjO TjO 

h T " H 0 
AH° 


spectroscopic constant for 
rotational stretching, cm - ' 

dissociation energy at 0° K, 
cal/mole 

Gibbs free energy for stand- 
ard state, cal/mole 

sensible free energy for 
standard state, cal/mole 

sum of sensible enthalpy at 
T° K and chemical energy 
at 0° K for standard state, 
cal/ mole 

sensible enthalpy for stand- 
ard state, cal/mole 

enthalpy change for forma- 
tion of substance from ele- 
ment in atomic gas state, 
cal/mole 




h 

Ah vap 

K 


Kf 


M 


N- 


n 


P 

Pi 

R 



enthalpy change for formation 
of substance from assigned 
reference element, cal/mole 

chemical energy at 0° K for 
standard state, cal/mole 

enthalpy, cal/g 

heat of vaporization to 1 gram 
of equilibrium vapor 

equilibrium constant for re- 
action of formation from 
element in atomic gas state 

equilibrium constant for re- 
action of formation from 
assigned reference element 

molecular weight based on 
chemical scale of natural 
oxygen, g/mole 

mole fraction of i Ln species, 

i-L 

moles of i Ln species per 
mole of mixture 

number of moles per gram of 
equilibrium mixture 

total vapor pressure, atm 

1 1 _ 

partial pressure of i in 
species, atm 

universal gas constant, 

1. 98726 cal/(mole)(°K) or 
82. 05971 (cc)(atm)/ 
(mole)(°K) 

entropy for standard state, 
cal/(mole)(°K) 
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s entropy, cal/(g)(°K) 

As vap entropy of vaporization to 

1 gram of equilibrium vapor 

T absolute temperature, °K 

v specific volume of vapor, cc/g 

1L 

x- weight fraction of 1 species, 

1 XU 

g of i species per g of 
mixture 

a vibration-rotation interaction 

e 

constant for diatomic mole- 
cule, cm' 1 

p density of vapor, g/1000 cc 

a standard deviation 


w 


e 


Gd x . y„ 
e e’ e J e 

Subscripts: 

c 

f 


g 

1 

2 


zero -order vibrational fre- 
quency for diatomic 
molecule, cm" 1 

anharmonicity constants for 
diatomic molecule, cm" 1 

condensed phase property 

formation from assigned 
reference element 

property of equilibrium 
vapor mixture 

monomer property 

dimer property 
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